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I. INTRODUCTION 


Despite the various astrophysical observations in support of its existence |l, 2], the nature of 
dark matter still remains an open question. Of the various candidates for dark matter, one of the 
most compelling is the Weakly Interacting Massive Particle (WIMP) 3-uf, with a mass ranging 
from a few GeV to 10 TeV. In fact, when the WIMP annihilation rate falls below the Hubble 
expansion rate, the chemical equilibrium between WIMPs and the primordial plasma is no longer 
maintained, and the number of WIMPs per comoving volume naturally fixes to the value required 
for explaining the present abundance of cold dark matter. Although chemical equilibrium at 
this stage is no longer maintained, kinetic equilibrium between dark matter and the plasma is still 
achieved through a high momentum exchange rate 8-3]. Eventually, when the Hubble rate equates 
the scattering process rate, WIMPs kinetically decouple from the plasma and flow with a given free- 
streaming velocity. This velocity sets the lowest value for the size of protohalos, which determines 
the subsequent evolution of primordial structures 17-23]. In particular, Bringmann [3] defined 


the temperature of the kinetic decoupling in the standard cosmological scenario, while Gelmini 
and Gondolo 21] defined T ^ in the Low-Temperature Reheating (LTR) cosmology following a 


dimensionality reasoning. 

In this paper, we present a full solution of the evolution equation governing the process of the 
kinetic decoupling, and we generalize the definition of the temperature of kinetic decoupling and 
the average kinetic energy of WIMPs in a generic non-standard cosmological model. This paper 
is organized as follows. In Sec. m we solve the evolution equation for the WIMP kinetic energy in 
a generic cosmological background, and we apply our results to the standard radiation-dominated 
scenario in Sec. uni to a power-law cosmology in Sec hyi and to a broken power law cosmology in 
Sec. El In particular, we discuss the Low Temperature Reheating (LTR) scenario in Sec. IV Bl and 
the kination scenario in Sec. IV Cl We summarize our results in Sec. El 


II. GENERAL SOLUTION OF THE TEMPERATURE EQUATION FOR DARK 

MATTER IN A THERMAL BATH 

The scattering process between plasma at temperature T and WIMPs of mass M x S> T can 
be described by a Brownian motion in momentum space. The momentum transfer A p and the 
average momentum of the dark matter particles p are related by 


P = V^Ve A P, 


(1) 
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where N e is the number of collisions required to change the momentum by p. Since p r^j \JM X T 
is much larger than the average momentum transfer A p ~ T, the number of collisions required 
to appreciably change the momentum of WIMP is N e = (jp/Ap) 2 ~ M x /T 3> 1. The momentum 
exchange rate T is suppressed with respect to the elastic collision rate T e i by a factor T/M x . 
Thermal decoupling of WIMPs occurs at a temperature Tkd approximative^ given by H(T^) ~ T, 
where H = H(T ) is the Hubble expansion rate at temperature T. A more precise way to define 


Tkj consists in using Boltzmann’s equation 


ion E, 


10-14 


to be treated as a perfect fluid. In particular, Ref. 


IQ], which does not require the dark matter 


14] discussed the Boltzmann equation under 


the generic assumptions of a heavy dark matter particle with M x T and with small momentum 
transfer per collision A p <C p, obtaining a Fokker-Planck equation for the dark matter particle 
occupation number f x = f x (p ), 


^ - H(T) p x • 
dt K ’ x 3p x 


= 7(T) ^'( p > / * (1±w+M ' T ls) 


( 2 ) 


The momentum relaxation rate q(T) is defined in terms of the scattering cross section between the 


• Q. 


Here 


WIMP and each relativistic species in the plasma at temperature T, see Eq. (2) in Ref. 
we are not interested in the specific form of 7 (T), since we are interested in a general solution to 
Eq. Q. We only assume that q(T) monotonically increases with T. 

Defining the WIMP kinetic temperature T x as 2/3 of the average kinetic energy of the dark 
matter particle, 




pj 

2 M v 


/x(Px) d3 Px> 


( 3 ) 


and using the Fokker-Planck Eq. ([2]) with the approximation 1 ± R ~ 1 and the Hubble rate 
H = a/a, it can be shown that T x satisfies the differential equation 


14, 


m 


da 


+ 2 T v = - 


27(0) 

H 


(T x -T). 


Eq. (j4|) is the central equation we consider in this paper. Refs. 


[ 12 ] and Q 


( 4 ) 


present analytic 


solutions to Eq. 
in Ref. 


in the standard cosmology and for 7 (T) proportional to a power of T, while 


. [[16] Eq. Q is solved numerically in the standard cosmology and in the presence of quark 
interactions. Here, we solve Eq. Q in terms of analytic expressions for a generic cosmological 
model. 

The scattering of the WIMP particles off the plasma is regulated by the function 

7(-n 


T(T) = 


H(T)' 


( 5 ) 
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in terms of which we rewrite Eq. (J4]) as 

dT 


da 


+ 2 [1 + T(T)] T x = 2T (T)T. 


( 6 ) 


The general solution for T x must satisfy the following boundary conditions. For y(T) H(T), 
WIMPs are tightly coupled to the plasma and the kinetic temperature approaches 


aT x = constant, 


( 7 ) 


from which, using Eq. Q, we obtain the behavior 


T — T v 


L x 


—>• 0, or T x ~ T, for T —> oo. 


x 2 T(T) 

In the opposite limit q(T) <C H(T), WIMPs decouple from the plasma, and Eq. (HD reduces to 


( 8 ) 


dT 


a 


x 


da 


+ 2 T x = 0, 


( 9 ) 


with solution 


a 2 T x = const. 


( 10 ) 


A. Analytic expression for the kinetic temperature T x 


We solve Eq. © by the method of undetermined coefficients, by first considering the homoge¬ 
neous equation associated with it, 

-,(hom) 


dT y 

a —^-+ 2 [1 + T(T)1 4 hom ) = 0, 

da A 


whose solution from the initial value of the scale factor a* to a is 

2 


( 11 ) 


j>(hom) _ ji(h°m) e ~ G (a,ai) 


( 12 ) 


Here, jl hom ) = T x loni \ai) is a constant temperature that will disappear from the full solution, and 


we have defined the function 


G(a,a') = 2 f T(a") 
Ja' 


da' 1 


(13) 


which satisfies the relation 


G(a, a') = G(a, a") + G(a ", a'), for a < a” < a'. 


(14) 
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Notice that, when expressed as a function of time, Eq. (1131) is simply given by 

G(t,t') = 2 [ 7 (t")dt". (15) 

Jf 

The particular solution to Eq. © is obtained by using the method of undetermined coefficients, 

Tj. part) (a) = Tf om ) (a) «(a), (16) 

where the function u(a) satisfies 

aTl hom \a) = 2T(a)T(a), (17) 

* da 

The solution to Eq. (HID is 

= A r g—G(a,a') a ' da '. ( 18 ) 

Ja.i 

The complete solution to Eq. (j 6 j) that satisfies the condition T x {ai) — T x i is the sum of the 
homogeneous and particular solutions, 

T x (a) = T X i (-Y e“ G(a ’ ai) + A e" G(a ’ a,) T(o') T(a') a' da'. (19) 

' a ' a Jai 

Eq. (|19D provides the value of the kinetic temperature T x as a function of the temperature of the 
Universe T for a generic cosmological model. Integration by parts of Eq. (I19|) using de~ G ^ a,a '^ = 
e ~G(a,a ’)2 T(a')da'/a' gives the alternative form 

T x (a) = T(a ) + (T xi - T t ) (A)' e~ G ^ - ^V G (“’“') ^7 A«') 


da. ( 20 ) 


Here Tj = T(oj) is the initial plasma temperature. 

Another alternative form of the solution is obtained by introducing the indefinite integral 

7 1 9 f a 7 A) da ' 

‘ {a) = 2 l (21 > 

which is defined apart from a constant that disappears from the solution. The integral in Eq. C3D 
can then be written as 


G(a, a') = s(o) — s{a'). 


Eq. (1191) in terms of the variable s = s(a) is rewritten as 


T x 


= T 


x* 




e s ~ s ' T(s') ds ', 


( 22 ) 


(23) 


where Sj = s(dj). 
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When the initial condition a* —>• 0, the exponential term in the homogeneous solution of Eq. (1191) 
drops to zero, while the product a* Tj remains constant because of the limit in Eq. (jSJ), yielding 


T x {a) = — f e G (“’ a ) Y (a')T(a') a'da', if a* —> 0. 

a Jo 


Equivalently from Eq. ([20]) . using T x i = T, at very small a^, 


T x(a) = T(a) -£e~^ ± [(^V) 


da', 


if cii -)• 0, T xi = Ti. 


(24) 


(25) 


Since a'T(a') —>• const for a' —> 0, the derivative in the integrand remains finite in the limit of tight 
coupling, while the exponential term drops to zero because of the limit in Eq. ([ 8 ]). Thus, T x ~ T 
in the tightly coupled limit. 


B. Temperature of kinetic decoupling 


The temperature of kinetic decoupling T kd expresses the temperature of the plasma at which 


the kinetic decoupling of WIMPs occurs. Here, we use the definition [l(j], 

7 (T kd ) = H(T kd ), 


(26) 


where i7(T kd ) is the Hubble expansion rate when WIMPs decouple kinetically from the primordial 
plasma. In the literature, different definitions of the temperature of kinetic decoupling can be 
found. In Bertschinger [12j], the definition of the temperature T^ d differs from our Eq. (1261) by a 


factor of two, 


'y(T&) = 2H(T&), 


- B N 


(27) 


while another definition of T kd is given by equating the rate of heat transfer equal to the Hubble 
expansion rate, 2y(T kd ) = i7(T kd ). 


Bringmann and Hofmann 


131 ] define the temperature of kinetic decoupling as 


rp2 

nnBH _ 1 

^ kd rp 

1 X 


(28) 


T-> 0 


Although the four definitions for T kd yield the same dependence on T x , up to a numerical constant, 
there are some conceptual differences between them. The definition in Eq. (|28|) depends in principle 
on the moment at which the temperature of KD is computed, and requires knowledge of the 
evolution of the universe at late times, far after kinetic decoupling. Instead, the definition we 
adopted in Eq. (126|) and Eq. (127]) depend on the properties of the WIMP-radiation coupling only, 
through 7 (T), and on the cosmology through H(T). As we show below, these expressions can be 
generalized to the case in which the cosmology is not in the form of a power-law model. 
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III. STANDARD COSMOLOGY 


In the standard radiation-dominated cosmology, it is common to parametrize the total energy 
density p(T ) and the total entropy density s(T ) in the universe with the so-called energy and 
entropy degrees of freedom g(T) and g s (T), respectively, defined so that 


7T 


p(T)=‘-g(T)T\ 


47 t 2 o 

s(T) = —g s (T)T 3 . 


The Friedman equation and entropy conservation then give the relations 

a 3 T 3 g s (T ) = const, 


H(T) = H iad (T) = T 2 J g(T), 


45 


(29) 


(30) 


where H iad (T) is the Hubble rate in the radiation-dominated cosmology. The scale factor depends 
on temperature as 


da dT { ldlng s \ 
~a ~ ~~T V + 3 din T ) 


It follows that 

G(T, T') = G(a, a') = 
and the WIMP kinetic temperature is 


,/ 45 , 

r' 7 cn f 

< 

to 

T T’*y/g(T") \ 


1 dinin' 
+ 3 d In T" 


(31) 


(32) 


T X = T 


, (T . _ T .) T 2 gs(T?' 3 ~G(T, Ti ) _ f T1 -G(T,T') T 2 g s (T) 2 / 3 ( 2 d\ug s [T’) \ 

{ XI ±l T 2 g s (Ti ) 2 / 3 J T T'2 5s (T')2/3 l, i + 3 din T' ) ‘ 


We have presented these formulas as they may 


dark matter numerical codes like DarkSUSY 


re useful when imp 


24] and micrOmegas 


(33) 

ementing these expressions in 


M- 


IV. POWER-LAW COSMOLOGICAL MODEL 

In this section we consider models in which the Hubble rate H, the scale factor o, and the 
momentum relaxation rate 7 have a power law dependence on the plasma temperature T. More 
specifically, we consider a dependence of the Hubble rate on temperature of the form 

H(T) = ^ , (34) 

where v is a positive constant, and T) and H{ are the temperature of the plasma and the expansion 
rate at the time at which we start considering the cosmological model. For the dependence of the 
scale factor on the temperature we write 


a a T = const. 


(35) 
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Equating Eqs. (IM1) and (|35|) . we obtain the relation 


H(a) = Hi - 
a 


(36) 


where dj is the scale factor at temperature 7). Notice that, in the radiation-dominated cosmology 
for which u = 2 and a = 1, the temperature of the plasma drops as T oc a -1 , while the WIMP 
temperature drops at a faster rate oc a -2 . For the momentum relaxation rate 7 (T) we assume 
a power-law function of the form 


l{T) = H 


I) 


4+71 


(37) 


where 7 * = 7 (Tj). In some models, the exponent n is related to the low relative velocity v of the 
forward WIMPs scattering amplitude _A4 forward off particles in the plasma, 


| M 


forward 1 2 


= const V 


Finally, in power-law models, Eq. ([5]) is given by 

4 +n—u 


T = — = T i 
H 


8 


— T 


a,' \ «(4+n-i/) 


(38) 


(39) 


where T ?; = y l /H l . 


A. Kinetic temperature 


Using the definition in Eq. (1211) in the power-law model, we find 


s a = 


2 T ? ; 


aj \ ot(A+n-v) 


a (4 + n — v) V a 


- 2 T, In ( — 

Qji 


, for 4 + n 7 ^ i/, 
for 4 + n = v. 


(40) 


In the first line, we have fixed the arbitrary constant in the definition of s(a) so that s(a ) is a power 
law. It is interesting to observe that for 4 + n/ y 


s = 


7 


a(4 + n — v) H 

Plugging Eqs. (l36l) . (l35|) . and (l22l) into Eq. (l23l) . computing the integrals, using the identity 


r(l + r, x) = r r(r, x) + x r e x , 


(41) 


(42) 


and defining 


A = 


2 — a 


a (4 + n — is)' 


(43) 
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we find 


T x = Ts x e s [r(l - A, s) + AT(—A,Sj)]. 


for 4 + n^i/, 


(44) 


and 


rr ( Q i \ 2 + 2T i 

T * = T ‘(j) + 


2T iT 


2 + 2T i - a 


1 - 


a,i \ 2+2T i-a 


(f) 


for 4 + n = v. 


(45) 


To the best of our knowledge, the expressions in Eqs. (|44H45I) have never been derived for the case 
of an arbitrary power-law model. In Appendix^! we show how to obtain the result in Eqs. (|44H45I) 
from solving the differential Eq. <[6|) directly. 

If the initial scale factor a* is taken so far back in time that the WIMPs are initially tightly 
coupled to the primordial plasma, then 7, 3> H, and .7 —» +00, and we obtain 

T x = TsVr(l — A,s). (46) 

Eq. (1461) is a generalization of the relation obtained in Ref. |l2] for any cosmological power-law 
model and for any value of the partial wave number n. 


B. Tkd in the power-law model 


Using the expressions in Eqs. (1341) and (1371) and the dehnition of the temperature of the plasma 
at which the WIMP decouples kinetically, Eq. (1261) . we obtain for 4-in^, 

Tkd = Ti (47) 

while in the case 4 + n = v we find the constraint 7j = Hi, for which no solution for T^d exists 
if this condition is not satisfied at all temperatures. In the literature, different definitions of the 
temperature of kinetic decoupling can be found. If we were to use the dehnition in Eq. (1271) . we 
would obtain 


T]f d = 2 l/(4 +n-v) Tkd> 


(48) 


Eq. (1531) below shows that the dehnition in Eq. (1481) is obtained by setting the variable s = 1 in 
the standard cosmology, for the case of a p- wave. 

Finally, the dehnition in Eq. (1281) applied to the power-law model discussed in this section gives 
T^ h = 0, unless in the particular case a = 1 to which radiation-dominated cosmology belongs to. 
For this reason, we suggest to modify the equivalent dehnition in Eq. (1281) as 


t vg _ T i 

J kd — TfT 


(49) 


T-> 0 
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Plugging Eq. (1461) for T x , valid for 4 + n / i/, in the definition in Eq. (1491) . gives 


t vg _ 

J kd — 


s~ x e~ s 


—X 


T, = 


F(l-A) 


r (i a, s) T ^ 0 

where Sj is a temperature-independent quantity defined by the relation in Eq. m as 

2 


Si = 


T 


(50) 


(51) 


a (4 4 -n — v) 

As anticipated in Sec. IIIBl the three results in Eqs. (1601) . (1451) . and (1501) have the same dependence 
on Tj, up to a numerical constant. However, the definition in Eq. (1491) depends in principle on the 
moment at which the temperature of KD is computed, unlike the definition we adopted in Eq. (1261) 
which only depends on the properties of the WIMP-radiation coupling through q(T) and on the 
cosmology through H(T), and can be generalized to the case in which the cosmology is not in the 
form of a power-law model. 


C. Radiation-dominated cosmology 


In the standard radiation-dominated model a = 1 and u = 2, Eq. (f43l) gives A = 1/(2 + n). 


Using the identity in Eq. (1421) . the limit for T x expressed in Eq. (1461) gives the same result as Eq. (4) 
in Ref. Q], 


T v = Ts 2 + n e SR T (s R ) = T 


2 + n 


1 -4-/1 

1 - e SR 4 +n T ( - 


2 + n 


2 + n 


sr 


(52) 


where the variable s defined in general through Eq. (1211) reduces in the radiation-dominated cos¬ 
mology to 

2 7 


sr = 


2 + n H 


(53) 


In addition, when we consider p -wave scattering (n = 2) in a radiation-dominated model, then 
sr = 7 / 2 H, and we obtain Eq. (12) in Ref. [jjj], 


T x = Ts « 4r | -,s R 


(54) 


D. Late time behavior 


When the plasma temperature is much smaller than T^, the late-time behavior of Eq. (|44|) gives 


T* = T iS 


x -i-i \ T 


T 


P(1 — A). 


(55) 
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In a cosmological model that approaches the radiation-dominated scenario where a = 1 and u = 2, 
Eq. (T55li reads 


T x = 


T 2 

~Ti 


2T,; \ ^ /l + n\ 
2 + n) \2 + n) 


We compare this result with the theoretical behavior 



rr ith 

X 


T 2 / 2 \ 2+^ /I + n\ 

2kd,std V^ + n / \2 + n/ 


(56) 


(57) 


where Tkd,std is the temperature of kinetic decoupling in the radiation-dominated cosmology, 


T+std = Ti 


H Tad (Ti)\ 2 + 
7 i 


(58) 


and H iad (T ) has been defined in Eq. (1301) . This latter equation can be stated in terms of the 
function T,; in Eq. (l60l) as 


T i _ rji 'Y' 

i — J- kd,std J- 


l 

2+n 


(59) 


This relation is also obtained by comparing the result in Eq. [56] with the theoretical Eq. (1571) . We 
rewrite Eq. (I59[) in terms of the temperature of kinetic decoupling Tkd by using the relation in 
Eq. (1601) in the form 

4+n —v 

, (60) 



as 


Tkd 


^^7u77,+2 ' 

^ kd,std 

rpis — 2 


4 -\-n — u 


= T t 


H rad (Ti) 

li 


4+n. — v 


(61) 


Eq. (1611) gives the temperature of the WIMP kinetic decoupling in a generic cosmological model, 
which might differ from the radiation-dominated scenario at the time of decoupling. Notice that, 
in the particular case in which the decoupling occurs in a radiation-dominated scenario (u = 2), 
Eq. (16T1) gives 


Tkd = Tkd,std- (62) 

In the following, we discuss the decoupling of WIMPs in a broken power law cosmological 
model, where a generic pre-BBN cosmology takes place before T t , after which standard radiation- 
dominated cosmology begins. 
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V. BROKEN-POWER-LAW COSMOLOGICAL MODEL 

A. Generic pre-BBN cosmology 


The relic density and velocity distribution of WIMPs depend on the characteristics of the 
Universe before Big Bang Nucleosynthesis (BBN), which occurred at temperature above Tbbn ~ 
1 MeV. Since this is an epoch from which we have no data, it is possible that the expansion rate 
of the Universe prior BBN differed from that of the radiation-dominated period, with the transit 


between t. 
above 


re two cosmological models occurring at a reheating temperature Trh which must lie 


26H3C] 


Trh > 4 MeV. 


(63) 


For this reason, in this section we generalize the power-law model discussed in the previous section 
by considering a generic scenario in which the Universe is dominated by some form of energy 
before cooling down to a reheating temperature Trh, after which standard radiation-dominated 
cosmology takes place. For this, we assume a dependence of the Hubble rate on temperature in 
the form 


H(T) = .Hrh 


(*)‘ 


for T > Trh, 


355 ) > for T < Trh, 


where v is a positive constant, Trh is the reheating temperature, and -£Trh = -H(Trh). 
scenario, the dependence of temperature on the scale factor is a power law of the form 


(64) 

In this 


T(a) = Trh 



for a < a RH , 
for a > a RH , 


(65) 


where orh is the scale factor when T = Trh- In the radiation-dominated cosmology the tempera¬ 
ture of the plasma drops asToc a -1 , while the WIMP temperature drops at a faster rate T x oc a~ 2 . 
Equating Eqs. ()64l) and (1651) . we obtain the relation for the Hubble rate in terms of the scale factor, 


H (a) = Hrh < 



for a < a RH , 
for a > a RH , 


( 66 ) 


Incidentally, a relation between Trh = 7(Trh)/-H(Trh) and T^ is obtained using the definition 
for the temperature of kinetic decoupling introduced in Eq. (I26p . in the form T(Tkd) = 1, and with 
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Eqs. ([571) and (1641) . 


Yrh 




4+n —v 
2+n 


for T kd > Tr H , 
for T k d < Trh- 


(67) 


Plugging Eqs. ([6511661) into Eq. (051) we find, for a < orh, 


?x=< 


T s x e s T (1 — A, s ), for 4 + n / u, 
for 4 + n = n, 


( 68 ) 


,2+2T rh — a ’ 

where s has been dehned in Eq. (141 p and A in Eq. (143 p . After reheating a > orr, we obtain 


m T 2 

T x — 7R 
Jrh 


sr 2 2+n 


e n S 


R 


El ^,s R )+C n 

2 + n 


where sr has been defined in Eq. (1531) . we set 

2 2 

SR = — Trh ’ s = -7 Trh ’ 

2 + n a (4 + n — v) 

and where the constant matching the two cosmologies at a = orh is 

Cn = S x se S - SR r (1 - A, s) - r , sr) , for 4 + n + /+ 


(69) 


(70) 


(71) 


and 


C - c Sr 2Trh 
" 2 + 2T rh — a V 2 + n 


1 H - 71 

— r [ --, sr ) , for 4 + n = n. 


(72) 


For high temperature T x ~ T, Eq. (1681) shows the limit discussed in Eq. (1461) . When the plasma 
temperature approaches T —> 0, Eq. (1551) gives 


T 


2 i 


Ty = —- S n 

x Trh r 


2+n 


r{ \r^ ]+Cn 


(73) 


Comparing the result in Eq. (f73l) with the theoretical expression in Eq. (|57|) gives 


i 

T i rp 'Y' 2+n 

RH — -L kd,std 1 RR 

This expression is similar to Eq. ([591) . with the additional inclusion of the matching constant C n 
which comes from the extra pre-BBN cosmology. 


1 + 


c„ 


r ( i±+ 

1 \ 2+n) 


(74) 
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B. Low temperature reheating cosmology 


In the Low Temperature Reheating (LTR) scenario 31-34], the energy density of a massive 
scalar field tfi dominates the evolution of the Universe in the pre-Big Bang Nucleosynthesis (BBN) 
epoch, after which standard radiation-dominated cosmology begins. The evolution of the Universe 
during the LTR stage is non-adiabatic: in fact, in the LTR model, the Universe is constantly 
reheated by the decay of the massive particle cj). This scenario is modeled by the Hubble rate in 
Eq. (1641) . with v = 4, and by the scale factor in Eq. (1661) . with a = 3/8. In the LTR scenario, 
Eq. (1741) for the reheating temperature as a function of Trh, with the constants in Eqs. (El) 
and d72l) . reads 


Tru — 

T kd ,std “ Irh 


1 + 


1 


r ( i±iA 

\2+n) 


s 


. 13 

S 3n e 
1 

Sr" e* R 


3n- 13 
3n ’ 


-r 


1 + n 

2 + n 


sr 


for n ^ 0, 

(75) 


and 


Erh _ -y-1/2 

T ~ 1 RH 

T k d,std 


1 + 


r® 


2Trh 


s R e Sli 


2 + 2Trr — a 


-r 


2’ 


sr 


for n = 0. 


(76) 


Fig. |T] shows the value of the temperature of kinetic decoupling T kc i as a function the reheating 
temperature Trh in the LTR cosmology, obtained from solving Eqs. ED and (1761) . Both temper¬ 
atures are given in units of Ekd,std- We considered the case of a dominating s-wave (blue), p -wave 
(blue), or d -wave scattering (red). Gelmini and Gondolo 21] recovered a relation between T kd and 


Trh in the case of the LTR cosmology and for p -wave scattering (see their Eq. (6) in Ref. 


T-tGG _ 
1 kd ~ 


L kd,std 

Trh ' 


for T kdjStd > Trh, 


21 ]) 

(77) 


where T kd std has been defined in Eq. 
Eq. ED with n = 2 gives 


Ekdst.d? for E kd]Std < Trh. 

Within the framework of the present paper, using 


E kd = < 


rp2 

1 kd,std 

Trh 

E kd)S td 


1 + 


r(|) 


1 + 


r ( 


C'2 


C 2 


for E kd > Trh, 
for E kd < Tr iH . 


(78) 


The main difference between Eqs. (17TD and (175D stands in the appearance of the constant C 2 . In 
Fig. [21 we show the relation between T kd and Trh for a dominating p -wave WIMP particle, using 
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FIG. 1. The temperature of kinetic decoupling Xkd as a function of the reheating temperature Trh in the 
LTR cosmology. Both temperatures are given in units of Tkd.std- We show results from Eqs. (ED and ED 
for s-wave (black), p -wave (blue), and d -wave scattering (red). 

the Gelmini-Gondolo model in Eq. (17TD (blue dashed line), and our result in Eq. (17HD (blue solid 
line). 


C. Kination cosmology 


In the kination cosmology scenario 


351 ], the expansion of the Universe before the standard 


radiation-dominated cosmology begins is driven by the kinetic energy of a scalar held. We model 
the kination scenario by setting the values for the pre-BBN cosmology a = 1 and v = 3 in Eqs. (1641) 
and dMD- With these values, the condition 4 + n ^ v is always satisfied for all n > 0, and the 
reheating temperature in Eq. (|74D in the kination cosmology is 


7r,h = Tkd.std { 1 + 


1 


r ( 

\ 2+n ) 


«!+" e 
1 

Sr + " e s R 


r 


n 


1 + n 


-r 


1 + n 

2 + n 


sr 


(79) 


Fig. [3] shows the value of the temperature of kinetic decoupling as a function the reheating 
temperature Trh hi the kination cosmology, obtained from solving Eq. (1791) . Both temperatures 
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ltr 2 ltr 1 io° 


TrH / Tkd.std 

FIG. 2. The temperature of kinetic decoupling Xkd as a function of the reheating temperature Trh, for a 
dominating p -wave mode in the pre-BBN LTR cosmology. Both temperatures are given in units of 7kd,std- 
Results are shown both using the Gclmini-Gondolo model in Eq. (1771) (blue dashed line), and the results 
derived in this paper from Eq. (1751) (blue solid line). 

are given in units of 7kd, s td- We considered the case of a dominating s-wave (black), p -wave (blue), 
or d -wave scattering (red). 


VI. SUMMARY 


In Eq. (124)) . we presented a general expression that gives the value of the WIMP kinetic temper¬ 
ature T x in terms of the temperature of the Universe T. This generic result has been specialized to 
the case of the standard cosmology in Sec. uni in view of possible future applications to numerical 
models. In addition, we have presented the expression for T x in the case of a power-law cosmology 
in Sec. IIV1 and for a broken power-law model in Sec. |Vj In Sec. IV B1 we have discussed the numeri¬ 


cal resu 
in Ref. 


ts in the case of the LTR cosmology, and we have compared our results with those obtained 


21] results, finding a discrepancy between our values of T^d as a function of Trh for low 


values of Trh- We have imputed this discrepancy to the appearance of the constant C n in Eq. (174|) . 
which come from the matching conditions between the LTR and the standard cosmologies. Results 
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FIG. 3. The temperature of kinetic decoupling Tkd as a function of the reheating temperature Xrh in the 
kination cosmology. Both temperatures are given in units of Tkd.std- We show results from Eq. (1751) for 
s-wave (black), p -wave (blue), and d -wave scattering (red). 

for the kination cosmology have been discussed in see Sec. IV Cl 
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Appendix A: Direct solution of Eq. (J6]) 


In the scenario described in Sec. m we solve the differential Eq. ([6]) by first writing T x as a 
function of the plasma temperature T, 

^ Tff + 2 + 7(T)) Tx = 27(T) T ' (A1) 

Here, dT/dt is obtained by taking the time derivative of T in Eq. (1351) and using Eq. (1341) . 

v+l 


— = -aHT = -aHiTi 


T 

Ti 


(A2) 







Using the result in Eq. (1A2I) and Eq. (l37l) . Eq. (1A1D reads 
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dT x 2 
dT a T 


1 + T,f s 


i \ 4 -\-n—u 


t y + 


2Tj fT\ A+n ~ u 


a \TiJ 


= 0. 


(A3) 


When 4 + ?i / we switch to the variable s defined in Eq. (1401) . expressed in terms of the 
temperature instead of the scale factor, 

4 +n—v 


s = s(T) = 


2 T ? ; 


a(4 + n-v) \Ti 


(A4) 


Rearranging Eq. (1A3I) . we obtain 


dT y 


ds \_a (4 + n — v) s 


+ 1 


TV + Ti 


a (4 + n — u) s 

2T 


4+n — v 


= 0, 


(A5) 


whose general solution is 

T x = r(j^j e s ~ Si +Ts x e s [r (1 - A, s) - T (1 - A, Si )] , (A6) 

Eq. (IA6D is identical to Eq. (144)) once the identity in Eq. (1421) is used. 

When 4 + n = v, Eq. (IA3I) reads 


dJ\ 

dT 


2 

a T 


Z 1 ; 

(1 + Tj) T x + —- = 0, 


(A7) 


and the solution is given by Eq. (1451) . 


T\ = T, 


<U\ 2+2 T 
a J 


+ 


2 Tj r 


2 + 2Tj - a 


1 - 


\ 


2+2T i-a 


(A8) 


To sum up, the solutions expressed in Eqs. (1A6I) and (1A81) . obtained by solving the differential 
Eq. (EH), are equivalent to Eq. (|44j) and (|45l) respectively, obtained by solving the integral Eq. (fl9|) . 
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